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Abstract 

We present a theory based on Green's function formalism to study magnetism 
in disordered Heisenberg systems with long range exchange integrals. Disor- 
dered Green's function are decoupled within Tyablicov scheme and solved with 
a CPA method. The CPA method is the extension of Blackmann-Ester ling- 
Beck approach to system with environmental disorder term which uses cumu- 
lant summation of the single-site non crossing diagrams. The crucial point is 
that we are able to treat simultaneously and self-consistently the RPA and 
CPA loops. It is shown that the summation of s-scattering contribution can 
always be performed analytically. While the p,d,f .. contributions are difficult 
to handle in the case of long-range coupling. To overcome this difficulty we 
propose and provide a test of a simplified treatment of these terms. In the 
case of 3D disordered nearest-neighbor Heisenberg system, a good agreement 
between the simplified treatment and the full calculation is achieved. Our 
theory allows in particular to calculate the Curie temperature, the spectral 
functions and the temperature dependence of the magnetization of each con- 
stituant as a function of concentration of impurity. Additionally it is shown 
that a virtual crystal treatment fails even at low impurity concentration. 
PACS numbers: 75.10.-b, 75.25. +z, 71.10.-w, 75.50. Cc 
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I. INTRODUCTION. 



The coherent potential approximation (CPA) is widely used to study the effect of disorder 
in crystals (for reviews see |l],|2|]). The CPA was initially developed independently by Soven 
|| and Taylor [fl to study systems with only diagonal disorder. Using a 2 x 2 formulation, a 
generalization to the presence of off-diagonal disorder was provided by Blackmann Esterling 
and Berck (BEB) In these approaches the main idea is to replace the system by an 

effective medium which is determined by the condition that the averaged T-matrix of a 
single impurity immersed in the effective medium is zero. An alternative approach is based 
on cumulant expansion |7|||. This latter method has the advantage that it can handle the 
environmental disorder term which is characteristic of the Goldstone's systems (phonons, 
magnons). The first proper treatment of the environmental disorder term, by using the 
cumulant expansion method is due to Lage and Stinchcombe who studied the diluted 
Ising problem (S=l/2). Later, using the 2x2 matrix method of BEB, the method was 
extended by Whitelaw |K| to the phonon problem. In their calculations the coupling and 
locator are fixed quantities and restricted to nearest neighbor exchange couplings. It is 
well known that magnetism in clean ferromagnetic systems can be tackled with Green's 
function formalism using Tyablicov decoupling procedure (RPA). This method goes beyond 
a simple mean field since it includes quantum fluctuations. Additionally, it fulfills the 
Goldstone and Mermin- Wagner theorems which is not the case of a mean field treatment. 
In the case of clean systems, combining first principle calculations to evaluate the exchange 
integrals and RPA method it was shown that one can provide satisfactory Curie temperature 



for Co and Fe Whilst, a simple mean field calculation largely overestimate the Curie 
temperature. It is our objective to provide in this paper a generalization of the RPA method 
to the disordered systems. We show that by combining in a self-consistent manner the RPA 
method and the CPA treatment of the disorder we are able to calculate Curie temperature, 
magnetization of the different constituants, spectral weights.... The CPA treatment is done 
in a similar way as done by Lage and Stinchcombe and by Whitelaw. However, due to 
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Tyablicov decoupling scheme for the disordered Green's functions, the locators and the 
effective exchange integrals are temperature dependent and have to be determined self- 
consistently for a given temperature. 

The paper is organized as follows. In the first section we derive after Tyablicov decoupling 
scheme the disordered binary alloy Green's function which includes diagonal, off-diagonal 
and environmental disorder. In section II, we perform the calculation of the averaged Green's 
functions for A (respectively B) atom. In section III, by generalizing Callen's formula we 
derive the equations for the magnetizations uia, trb and for the Curie Temperature. In 
section IV, we propose an alternative simplified treatment of the p,d,.. scattering contribution 
to the self-energy to the case of system with long-range exchange coupling. Finally in section 
V we present some numerical results and proceed to a test of our approximation of the self- 
energy contribution of the higher scattering terms. 

II. DISORDERED GREEN'S FUNCTION AND RPA DECOUPLING SCHEME. 

We study the magnetism in a binary alloy Ai- C B C , A and B can be either magnetic ions 
or paramagnetic. We denote their spin respectively Sa and Sb- The total Hamiltonian 

reads, 

H = £ - JyS* -Sj-Y: Di(S*) 2 -BY, 9l*{Sf) (1) 

ij i i 

where the Jy and Di are random variables: Jy = J^ y with the probability P^Pf where P t x 
is the probability that the site i is occupied by a A-atom. Similarly Di = D\ with probaility 
Pj A . The exchange integrals are assumed to be long range, our study is not restricted to 
the nearest neighbor Heisenberg model. The second term which describes anisotropy is only 
relevant in the case of 2D systems to get a non zero Curie temperature T c (Mermin- Wagner 
theorem). However in the case of 3D systems the contribution of this term can be neglected. 
We also include the effect of an external magnetic field. 
Let us consider the following retarded Green's function, 
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G±-(t) = -i9(t)([Snt),Sr(0)]) 
where (..) denotes the statistical average at temperature T, 

(6) = ^Tr(e-^O) 

where Z = Tr(e _/3H ). 

G^~(t)'s Fourier transform in Energy space is, 

/+oo 
-oo 

Its equation of motion reads, 

uG±-(u) = 2m i 5 ij + < [Sf, if]; Sr > 

where m; =< Sf >, or m; = (resp. tub) if i — A (resp. i = B). 

After expanding the second term on the right side of the equality we obtain, 

(u - gniB)Gtf{uj) = 2rrii5ij 

y t J-ii Si st st Si ; Sj 3> 
I 

+A < SfSt + S+Sf > 



(2) 



(3) 



(4) 



(5) 



(6) 



The next step consists in decoupling the higher order Green's function. For the second 
term we use the standard Tyablicov decoupling |12[ (equivalent to RPA). The last term due 
to anisotropy is somehow more complicate since on-site correlation are involved. Following 
the approach discussed in Ref. [[Uj we adopt for this term the Anderson-Callen decoupling 
scheme |14) |: 



A < S?S? + S*St >= 2Da imi 



(7) 



where, 



7< = 1 



2S 2 



{s i {s l + i)-<{siy> 



(8) 



After simplification we find, 



Gtf = giSij + 9i J2 ®uGt ~ w(E *») G i~ ( 9 ) 
i i 

where 4>u = — 1/2J# and = —l/2Ju^ and ^ denotes the locator: g t = g\ (resp. g B ) 
if i = A (resp. i = B). 

^ = E-^B/Zn-D^ (10) 
where A = A or B. For convenience, we have also introduced the reduced variable E = 
m denotes the averaged magnetization: m = Y,\ c \ m \- The term which is proportional to e 
comes from the environmental disorder term. This term is crucial to recover the Goldstone 
mode and requires to be treated very carefully. We have introduced the coefficient e which 
is in principle equal to 1, in order to follow the influence of the environmental disorder term 
during the calculations. Note also that this term appears because of RPA decoupling. If 
e = Eq. H is analogous to the propagator of an electron in a disordered medium with 
on-site potential and random long-range hopping terms tu = <&u [off- diagonal disorder). In 
this case the problem can be solved just within the BEB formalism. However, one should 
stress that the BEB formalism does not apply when the environmental term is present. Note 
also that in our model the locator g®, ^ u and ji are all temperature dependent, thus CPA 
and RPA loops have to be treated simultaneously in a self-consistent manner. It is also 
interesting to note that ^ in the case where the sites i and 1 are occupied by different 
type of atoms. 

III. CUMULANT EXPANSION METHOD FOR THE AVERAGED GREEN'S 

FUNCTIONS. 



As it is done in Ref. [HJ, the basic idea is to write Eq.|9| as a locator expansion in BEB 



manner ||. We define the random variable pf p$ — 1 if A is at site i or pi = if i is occupied 
by a B ion. Therefore the locator reads, 

g l =p l g° A + (l-p l )g° B =gf + gf (11) 
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and, 



PiJ$ A Pi+Mr(i-pi)+ 



i-Pdjrpi+^-PiWi^-pi 



(12) 



similarly, 



PiJi A Pi+PiJu a '\i-pi) + 



ABA, 



i- Pi )4 B > 2 Pi+(i-Pi)jr(i-pi 



BB i 



(13) 



where J„ 



AB,1 

il 



^J t AB and 4 B ' 2 



m A jAB 



m B 



a 



The Green's function are expressed in term of a 2x2 matrix and one gets for the equation 
of motion, 
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(14) 



tBB 



V jBB 
1^1 J il 



We have defined the variables J AB ' 1 = J2i Ju ' an d J L — i^i o u 
The aim is to expand this expression into a product of the p factors, which can then 
be averaged over disorder by expanding into cumulants. For that purpose we separate out 
the factors and introduce a new variable pi by pi = pi + c (where ca = c) . The idea is to 
separate out the virtual crystal part. 



(15) 



There is still the environmental term which is more difficult to handle. As it was done 
by Lage and Stinchcombe ||, by converting into k-space the calculations become easier to 
perform. 

We define the Fourier transform by, 
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G kk / = exp(ik • Ti) exp(-z'k' • r^G; 



'.i 



After some manipulation one gets 

/ 

kk' — Pk-k' + 



C 



\ 



°~ 1 / 

where the 2x2 matrix Vk q is defined by 



^k-k' + ^*k° X! Pk-qVkqG q k' 



V 



kq 



( T AA - c ( T AA - T AB '^ 
■ J q fc Wk-q J k-q J 



AB 



\ 



jAB 
J 1 



J. 



jBB i jBB j A B,2\ 
J q H J k-q J k-q ) J 



and the virtual-crystal Green's function G k c , 



[G k c ] = M - cMi 



where the matrices M and Mi are, 



Mn 



(9. 



A)' 1 







\ / eJ AB,l 



V 



IB) 



+ 



jAB jBB , jBB 
J k J k tJ 



and, 



( J AA - e(J AA - J^' 1 ) 

JBB _ t{J BB _ jAB,2 

The equation (|i~7f) can be expanded into 2 sub-series. 



Mi 



J k 



V 



tAB 
J k 



*^kk' - ^kk' + ^kk' 



(16) 



(17) 



(18) 



(19) 



(20) 



(21) 



(22) 



where the sub-series are respectively, 



G k k' - GJfpk-k' + — Y Gk C V kq G q C p k -qPq-k' + •••• 



N 



and, 



G k k' - ^G k c 5 k -k' + G k c Vkk'G k ^p k -k' + Gk 6 V kq Gq C V q k'G k ^pk-qPq-k' + 



(23) 



c 
c- 1 



(24) 
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The averaged Green's function is obtained by averaging over products of p by expanding 
into cumulants Pj(c). For instance, 



(Pk 1 p k2 ) = ^^(k 1 + k 2 ] 



N 



(pidPkaPka) 



A(c) 
N 2 



<J(ki + k 2 + k 3 ) 



(25) 
(26) 



and, 



(P kl Pk 2 Pk3Pk 4 ) 



P 4 (c^ 



5(ki + k 2 + k 3 ) 



[5(kj + k 2 )<5(k 3 + k 4 ) + 



<J(ki + k 3 )<5(k 2 + k 4 ) + 5(kx + k 4 )5(k 2 + k 3 )] 



(27) 



The cumulants are systematically obtained the generating function, 



g(x, c) = ln(l - c + ce x ) = ]T Pi{c)%r 



(28) 



i=l 



From this equation one gets Pi(c) = c, P2{c) — c(l — c), ^(c) = c(l — c)(l — 2c).... 

In order to get a closed form for the series we have to make the usual CPA approximation 
which consists in keeping only the diagrams with no crossings of external lines. As it is was 
shown by Yonezawa et al. J7||§, the self-consistency requires a modification of the semi- 
invariants to be attributed to each vertex. In other words it means that the cumulants Pi{c) 
have to be replaced by a new set of coefficients Qi(c) which satisfies the relation, 



Qi(c) + Q 2 (c)x + Q 3 (c)x 2 = a c (x) 



1 —x(l - cr c (x)) 



(29) 



where the modified cumulants are, 



m=l 



(i + m-2)! 
m\(i — m)\(m — I) 



(30) 



In the single site approximation, after averaging, one gets for the averaged 2x2 Green's 
function matrix, 



'kk' 



r k"k-k' 







/ 

c 


o ^ 




k' — Gk 








+ A k 








c-lj 





(31) 
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where, 

(32) 

£ k denotes the self-energy, it is given by, 

S k = Q 2 ^ E V kq G q V qk + Q 3 ^ E V kq G q V qt G t V tk + ... (33) 

q iV q,t 

and, 

A k = Q 2 i E V kqG q + Q 3 ^ E V kq G q V qt G t + ... (34) 
q q,t 

The term A k which is very similar to the self-energy is called end correction |1. Note 
that, inside the CPA loop, Eq.(|3"2"D and Eq. fl3"3]) are the only 2 equations which have to be 
solved self-consistently. To summarize, in Fig. [TJ we show a diagrammatic representation of 
the previous set of equations. 



G, 



[GIT 1 



A. Evaluation of A k . 



It is convenient for the calculations to start by defining, 



T«(q) = — E ex p( i( i r ?) ( 35 ) 

Zi 



The sum, t\ runs over the i-th type of neighbors of the i-th shell from a given site and 
Zi is the total number of neighbors in the shell. Note that, from now J2i will correspond to 
a summation over the different shells. With this definition it follows immediately, 

^(q)=E^^(q) (36) 

i 

We get similar expression for J BB (q) and ^^(q)... 

It is convenient to decompose the matrix Vk q into two terms, 

V kq = Vg + Vg (37) 

where 



Vg = E Vkq' = E [A* - eD ^(k)] 7< (q) 



r(l),i 



(38) 



and, 



Vg = E V kq" = [ 7< (k) 7 i(q) - 7,(k - q)] 



r(2),« 



(39) 



Aj and Dj are the following 2x2 matrices: 

/ 

A i = 



A.\ I A3 \ 



1 AB JP B 



(40) 



( jAA _ jAB. \ 



D, 



V 



o J BB - jf B ' 2 



(41) 



By using the following very useful property |T5|] : if f(r) is a function which is equivaluated 
at each site r,- of £7,- then, 



1 £ 7i (k - q)/(q) = 7(k)^E 7i(q)/(q) 

q q 



(42) 



By using Eq.(^), we find significant simplifications in the calculations. Indeed all the 
terms of the serie involving at least one factor reduces to zero. Thus the end correction 
term does not explicitely depend on the environmental disorder term. 

After calculation we finally get, 



A k = Evg,' 1 g 2 I + Q 3 M + Q 4 M 2 + .... 



(43) 



Like V^ 1 -*' 1 , F 3 is a 2 x 2 matrix, and M a N s x N s matrix, where each matrix element 
Mjj is a 2 x 2 matrix. N s denotes the number of considered shells. V^ 1 ^* is given in Eq. 



and F- 7 and My are defined by 



Fi = T7E7,(q)G q 

q 



(44) 



and, 



1 



M ^ = ^E7,(q)G q VS 
q 



(45) 
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The sum in Eq. [43] is obtained after diagonalization of the 2N S x 2N S matrix M = 

Q 2 I + Q 3 M + Q 4 M 2 + .... = p- 1 [( ( r c (M dia9 ) - Q^M^P (46) 

The function er c was previously defined in Eq. (|29"P, and [cr c (M. diag )}ij = cr c (Aj)<5jj where Aj 
are the eigenvalues of M. Hence, we get for the end correction 

A k = E vS J [p- 1 [(a c (M dM9 ) - QflM&jPhFi (47) 

Let us now proceed further and evaluate the self-energy £ k . 



B. Evaluation of £ k . 

Using the remarks made in the previous section, we find that the self-energy can be 
written, 



(i) _l 

k 



(48) 



where S k (resp.S k ; ) is obtained by replacing V kq by V k ^ q (resp.V k ^ q ). Indeed we 
find that each term of the serie containing both and V^ 2 ^ reduces to zero. After 



< 2 h ; 



(i) 



,(2) 



simplifications we obtain for S k 



S (1) (k) = v Vo toil + Q 2 M + Q 3 M 2 



1, j 



..r j (k) 



(49) 



where r j (k) = 7j(k) 



As previously done for the end correction, using the function a c (z) defined in Eq. 
obtain immediately, 



wc 



Q X I + Q 2 M + Q 3 M 2 



p->(M dw ,)]P 



(50) 



Note that, we have included in the sum the first order term depending on c (Qi) which 



comes from the virtual crystal Green's function G 



vc 

g ' 
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In general, the evaluation of the second term I]( 2 )(k) is much more complicated. One 
can get an analytical form only for simple cases. For example if the exchange integrals are 
restricted to only nearest neighbor, the complete summation of the serie can be performed 
by using the space group symmetry of the lattice fllBUj. In the case of nearest neighbor 
Heisenberg system one gets, 



4 2) (£) = C p (l - 7 (2k)) + C d (l + 7 (2k) - 2 7 (k) 1 



(51) 



where, 



: (Q l I + Q 2 M P;d + Q 3 M 2 Ptd +...)r> 1 



(52) 



C Pt d are evaluated in the same way it was done for Y,^\e) and A.^(E). The matrices D 1; 
M p and M d are respectively, 

( t A A T 4R1 ^ \ 



Di 



jAA _ jABA o 

JBB _ JAB,2 



(53) 



-DlGr 

6 ? 



(54) 



(55) 



where G p = ± £ q (l - 7(2q))G(q) and G d = ± £ q (l + j(2q) - 2 7 (q) 2 )G(q). 
Note that the virtual crystal approximation for ~S^\e) consists in taking in Eq. ( 52 ) 
the first term only. Then it follows immediately, 



rtvuA _ n 



VGA 



— D x 

2 



(56) 



which substituted in Eq.(^TJ) leads to, 



^' VCA (E) = -ecB 1 (l- 1 (kf 



(57) 



Note that £ 



(2),VCA 



is energy independent. It is also important to stress that at the lowest 



order the self-consistency for X^ 2 ) is not required. 
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Most of the ferromagnetic materials are of itinerant type, which means that the exchange 
integrals between different localized magnetic ions are long range and driven by the polar- 
ization of the conduction electrons gas as it is for the RKKY mechanism [D] . Analytically, 
the generalization of the previous calculations to the more interesting case where Jij are 
long ranged is not an easy task. However by truncating the serie, the summation can be 
performed numerically. It is important to note that S^(k) is (i) proportional to e which 
means that it originates only from the environmental disorder term and (ii) each term of the 
serie vanishes in the long wave length limit S^ 2 -*(k = 0) = 0. This implies that even after 
truncation of the serie at any order, the Goldstone theorem remains fulfilled. Thus the long 
wave length magnons are always treated properly. Furthermore, since S^(k) corresponds 
to higher order scattering terms (p,d,f,...) it is natural to expect that these terms should not 
affect the Curie temperature in a dramatic way. In other words we expect that a truncation 
of S( 2 )(k) serie to the first few term should already provide a good approximation of Curie 
temperature compared to the one one would get with the complete series. However, it is 
crucial to consider at least the lowest order term (the virtual crystal contribution) other- 
wise even in the clean limit one would not recover the correct result and the Goldstone 's 
theorem would be violated . If we consider the lower approximation 5^2) ~ ~Sy CA , we get 
the expected results in the limit c = and c = 1. It is not a priori clear whether such an 
approximation of S^ 2 ^(k) to the lowest order provides satisfying results for the Curie tem- 
perature at moderate impurity concentration. Such an approximation will be tested later 
on. 

To conclude this section the complete averaged 2x2 Green's function is obtained after 
solving self-consistently the set of equations Eqs. (|32|) and ( p3|) within the CPA loop and 
then using Eqs. ^1] and |33| to get and Gk- However, as was already mentioned in the 
introduction, the problem is not solved until we are able to calculate the locators and the 
exchanged integrals ^>u which depend on the averaged magnetization m\. The determination 
of m\ has to be done self-consistently in an additional external loop (RPA loop). 
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IV. MAGNETIZATION AND CURIE TEMPERATURE. 



We assume that the averaged 2x2 Green's function matrix G(k, E) is calculated ac- 
cording to the previous section within the CPA loop. We show how from G\(\i, E), X = A 
or B we can get the missing self-consistent equations (RPA loop) to get the temperature 
dependent locator g\ and the exchange integrals ^u. This will allow us to calculate the 
element-resolved magnetizations m\ =< S% > as function of temperature and the Curie 
temperature. It was shown by Callen, in the case of a clean system (pure A or B) that the 



magnetization can be expressed in the following way |15] , 



(Si - + * A )' i5 ->+ 1 + (Sj, + 1 + *i)*f * +1 ,„. 
mA = ( i + g >P> +i < 58 » 



where $a = ^Sq'^q) and $a(q) is defined as, 



where, 



I yE /j E % E } i (59) 



A x (d,E) = — ImG+-(q,B) (60) 

7T 



is the spectral function. 

Note also that the Callen's approach to get the magnetization allows to derive a lot of 
local spin-spin correlation, they are only expressed as a function of $a- For instance, 

((S z x ) 2 ) = S(S + l)-m x (l + 2$ x ) (61) 

which is needed to determine the anisotropy parameters 7a given in Eq. (|8]). 

In the case of clean systems, the normalized spectral function A\(q, E) is given by 

A x (q,E)=5(E-E(c l )) (62) 

E(q) = w(q)/2m and c<j(q) denotes the magnon dispersion. 

In the case of a binary (or multi-component) alloy this formula can be generalized in the 
following way, 
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where c\ is the concentration of A-ion and we have for convenience introduced a In- 
dependent reduced variable x\ = 

Note that in the presence of impurities the spectral function is not anymore a 5 function, 
but because of the finite imaginary part of the self-energy will consists of peaks of finite 
width with more or less a Lorentzian shape. In the case of binary alloy we expect for a given 
q, 2 peaks, more generally n peaks for an n-component alloy 

For a given temperature the complete self-consistency is obtained by (i) providing good 
starting values for m\ then (ii) performing the CPA loop which provide G(k, E) and finally 
(ii) going into the RPA loop by using Eq. (|58D , fl6~l|) and (^) one gets the new values of rri\ 
and ((Sa) 2 ) which are re-injected in the locators g®, the exchange integrals and 7a- 

Let us now show how to get the Curie temperature of a disordered Heisenberg binary 
alloy. We start by expanding Eq. (KO) in the limit T —>■ T c (i.e, m\ — > 0). We immediately 



get, 



where, 



q 

After expanding Eq. |58] CIS db function of one obtains, 



kT 

$A « f^A (64) 
2m 



i-EP^y (65) 



S A (S A + l)2m 1 

m > = ^^w c r x (66) 

Since the averaged magnetization m is defined by, m = J2\ c \ m \, combining the two 
previous equations one finds for the Curie- Temperature, 

7 rp 2 ^ S\(S\ + 1) 

k B% =^2^ c a p (67) 

6 A -\ 

Eq. (|67[) is the RPA generalization of the Curie Temperature to a multi-component 
disordered alloy. The previous equation provides a direct measure of the weight w\ = 
^^-[cA^fr^] of each A-element to the Curie Temperature. 
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V. NUMERICAL RESULTS. 



In this section we provide an illustration of the RPA-CPA theory and a test for the 
approximation suggested above for the higher order scattering contribution of the self-energy. 
For simplicity, we consider the case of a 3D disordered binary alloy on a simple cubic lattice. 
Additionally we restrict the exchange integrals to nearest neighbor only which allows us to 
test the validity of the approximation scheme suggested in Sec III Bfor estmating E 2 . For 
further simplifications of the calculations we consider the cas of a zero external field and 
neglect the anisotropy term which is reasonable for a 3D systems. 

In Fig. [|, we have plotted the Curie Temperature as a function of c obtained with the 
full CPA treatment, the E^ part of the self-energy is calculated exactly (full summation of 
the serie). Note that pure A (resp. B) corresponds to c = 1 (resp. c = 0). Depending on 
the chosen set of parameters T c shows ( i) a minimum (J abS 'aS b — min(J t aaS\, JbbS"b) , (ii) 
a maximum (JabSaSb > max(J 'aaS\, JbbS%) or (iii) is monotonic (min(J aaS\, JbbS%) < 
JabSaSb < max^AAS*!' JbbS%). These three different cases are shown in the figure. 

As already mentioned in section III, it is difficult to perform the full summation of E^ 2 ) for 
the case of long-range exchange integrals which is the case of many realistic and interesting 
systems, for example permalloy. As it was discussed previously the simplest approximation 
consists in keeping only the lowest order term of the serie ( virtual crystal approximation). 
In the case of nearest neighbor Heisenberg system E^ 2 ) and YP>> VOA are respectively given 
in Eqs (|51~D and (57) . In Fig. |] we have plotted the Curie Temperature calculated with (i) 



a full CPA treatment, (ii) the one performed with the approximation discussed previously 
$]( 2 ) = "Eyc A and (iii) the one obtained with virtual crystal approximation. In the case (iii), 
the averaged Green's function is, 



G k = GJf 



<c \ 



V 0c-1, 



(68) 



since in VCA A k = 0. 

The comparison between the full CPA and the virtual crystal approximation shows that 
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the Curie temperature differs significantly. Even, very close to the clean limit the VCA 
appears to be inappropriate, for instance for c = 0.1, we observe that T^ CA is about 35% 
larger than the full CPA calculated one. Note that the disagreement is even more pronounced 
in the vicinity of c = than c = 1. This can be understood in the following way: since 
Jab = 3Jbb = 1-5 Jaa an d Sa = Sb a substitutition of a B site by a A site (close to 
c = 0) introduces a change of energy (with respect to the pure case) 2 times larger than a 
substitution of a B site by a A site near c = 1. As discussed previously it is interesting to 
compare the Curie temperature where the VCA is only done on S*- 2 -* (T 2,VXM ) . We observe 
a good agreement between the full CPA calculated T c and T^ ,VCA , in the whole range of 
concentration, the agreement is even excellent for c > 0.6. A comparison between T^ CA and 
rp2,vcA - n ^ e vicinity of c = and c = 1 shows that the reason why the VCA approximation 
breaks down is essentially because of the crude approximation of the s-part of the scattering. 
Thus this figure validates a simple treatment of S^. It is also expected that including only 
few additional terms of the serie will lead to an excellent agreement in the whole range of 
concentration. Additionally, the approximation S( 2 ) « ^( 2 >> VCA will get better in the case 
of long range exchange integrals. 

In Fig. [| we show the temperature dependence of the element-resolved magnetizations. 
In order to demonstrate the versatility of our approach, we have chosen a set of param- 
eters which mimics a ferrimagnetic behavior with compensation point. Additionally, the 
parameters are such that T A 3> T C B . While the temperature dependence of follows a 
standard behavior, mg(T) start to strongly decrease even at low temperature. For example 
at T « 2.5, ttla has reduced by less that 20% while m# = 0.5ms(0). As a result of our 
choice of the parameters we see that the averaged magnetization m av = Icattia + c^m^l is 
non monotonic and vanishes for an intermediate temperature value (compensation point). 
It is found that the function ^-(T) decreases monotonically with temperature. As a result 
and since at T = 0, 2p = thus if ^ < ^ then m av will not have a compensation point 
. However, the condition that I s > — is not sufficient to get one, it also required that 
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mA \ o — Cb 

In Fig. |5| we now show the magnon spectral density (MSD) p\(E) = }^ as a function 
of E. We consider 3 different cases: almost clean A and B ( (a) and (c)), and the intermediate 
situation ca = cb = 0.5. In both, Fig. [5]a and Fig. |5]c we observe that the MSD is very 
similar to the clean case. This is clearer in case (c) than (a), it is easy to understand that 
when doping A with B the difference in energy with the undoped case is only of order 10% 
(Jaa(Sa) 2 = 0.8 and JabSaSb = 0.9) whilst doping B with A the change is more drastic 
(about 100%). To get a similar MSD to Fig.|]c for a weakly doped B sample, one should 
take c » 0.005. 

In Fig. |6] we show the spectral function S\(q,E) as afunction of energy for different 
values of the momentum q. This quantity is more interesting that the integrated MSD since 
it provides direct information about the elementary excitation dispersions and their spectral 
weight. Additionally it is directly related to inelastic neutron scattering measurements. Let 
us now briefly discuss Fig. |6|. At precisely q = momentum, in both S\=a,b, we observe 2 
peaks structure (i) a well defined peak |Tj| at E = as expected since our theory fulfills the 
Goldstone theorem and (ii) a very broad one at intermediate energy E w 1. For intermediate 
values of the momentum, it is difficult to separate the peak and one get a single broad peak. 
We see clearly that the peaks are crossing each other at q ~ |(1,1,1). Note that due 
to (i) the different spectral weight of the peaks and to the closeness of their location, the 
single peak-structure which is observed at q = ~ is located at different energy for A and B. 
From this figure we see also that the dispersion of the second peak is almost flat £"2(9) ~ 1, 
while the Goldstone mode -E'i(q) [0 goes from E = to E max m 2 when moving in the 
(1, 1, Indirection. 

VI. CONCLUSION. 

In conclusion, we have presented in this paper a theory based on Green's function for- 
malism to study magnetism in disordered Heisenberg systems with long range exchange 
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integrals. The disordered Green's function are decoupled within Tyablicov procedure and 
the disorder (diagonal, off-diagonal and environmental) is treated with a 2 x 2 modified 
cumulant CPA approach. The crucial point is that we are able to treat simultaneously and 
self-consistently the RPA and CPA loops. Our theory allows in particular to calculate Curie 
temperature, spectral functions and temperature dependence of the magnetization for each 
element as a function of concentration of impurity. Additionally we have proposed a sim- 
plified treatment of the p,d,f .. contribution of the self-energy which is difficult to handle in 
the case of long range exchange integrals. The approximation was tested successfully on 3D 
disordered nearest-neighbor Heisenberg systems. Combined with first principle calculations 
which can provide the exchange integrals this method appears to be very promising to study 
magnetism in disordered systems. 
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FIGURES 

FIG. 1. Diagrammatic representation of the averaged Green's function calculated within 
the CPA loop. G is the total averaged Green's function, is the self-energy and the 
end-correction. 

FIG. 2. Curie temperature T c for disordered nearest neighbor Heisenberg ferromagnet as a 
function of the impurity concentration c (A). The parameters are Sa = 2, Sb = 3, Jaa = —0.2 
and Jbb = —0.15. We have chosen 3 different values for Jab- 

FIG. 3. Comparison between the Curie temperature calculated as function of the impurity 
concentration for a nearest neighbor Heisenberg Ferromagnet. (a) the full CPA calculation, (b) ap- 
proximation for E^ 2 ) = ^vca an d ( c ) the virtual crystal calculation. The chosen set of parameters 
are written in the figure. 

FIG. 4. Magnetizations tua, —tub and averaged one c av = \cAfnA + Q>m.B| as a function 
of temperature. The spins are Sa = 1 and Sb = 3, the exchange couplings are Jaa = —1-2, 
Jbb = —0.10 and an anti- ferromagnetic coupling between A and B is taken Jab = 0.15. The 
concentration of A-atoms is ca = 0.70. 

FIG. 5. Density of state p\(E) = as a function of E. The continuous line corresponds 

to A = A and the dashed line to A = B, for 3 different concentration of A: c = 0.05, 0.5 and 0.95. 
The parameters are Sa = 2, Sb = 3, Jaa = —0.2, Jbb = —0.05, Jab = —0.15 and T T c . 

FIG. 6. Spectral function S\(q,E) = -±ImG x (q,E) as a function of E for different momentum 
q where q = q(l, 1, 1). The continuous line corresponds to A = A and the dashed line to A = B. The 
spins are Sa = 2 and Sb = 3, the exchange couplings are Jaa = —0.2, Jbb = —0.10, Jab = —0.15 
and ca = 0.50. We have taken T « T c . For clarity of the picture a small imaginary part r] = 0.1 
have been added. 
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